High-temperature series expansion of the Helmholtz free energy of the quantum 

spin-S XYZ chain 
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We consider the XY Z chain model of arbitrary spin S in the high temperature region, with 
external magnetic field and single-ion anisotropy term. Our high-temperature expansion of the 
Helmholtz free energy is analytic in the parameters of the model for S, which may range from 1/2 
to the classical limit of infinite spin (S — > oo). Our expansion is carried out up to order (J/3) 5 . 
Our results agree with numerical results of the specific heat per site for S = 1/2, obtained by 
the Bethe ansatz, with h = and D = 0. Finally, we show that the magnetic susceptibility and 
magnetization of the quantum model can be well approximated by their classical analog in this 
region of temperature. 



The one-dimensional Heisenberg model describing the 
exchange interaction of spins is intensively applied to 
study the magnetic properties of quantum spin chains. 
The XXZ spin- 1/2 chain with external magnetic field 
along the z direction is exactly integrable, consequently 
its thermodynamics has been investigated since very 
early works P, 0- It has also been verified the ex- 
istence of higher-spin quasi one-dimensional magnetic 
systems such as CsVCl 3 and CsVBr 3 (S = 3/2)0, 0, 
(CioHgN^MnCls (S = 2)@ and (CD 3 ) 4 NMnCl 3 (S = 
5/2) U0- AH these antiferromagnetic structures exhibit 
nearly ideal one-dimensional behavior over a consider- 
able range of temperature. Motivated by these features, 
magnetic and thermodynamical properties for higher- 
spin chains were first investigated numerically 0- More 
recently we obtained in Ref. |i| the /3-expansion of the 
Helmholtz free energy (HFE) of the XXZ model, for ar- 
bitrary spin-S, up to order (J/3) 6 , where /3 = 1/kT (k 
is the Boltzmann constant and T is the absolute tem- 
perature), in the presence of an external magnetic field 
and a single-ion anisotropy term. We showed that in the 
high temperature region the magnetization and the mag- 
netic susceptibility, even for low values of spin such as 
S = 3/2, can be properly approximated by the classical 
results within a percental error of less than 4%. 

Quasi-one-dimensional magnetic systems have 
been synthesized experimentally such as the 
(C 6 H 11 NH 3 )CuBr 3 and the (C 6 H 11 NH 3 )CuCl 3 , known 
as CHAB and CHAC, respectively |lCj, both exhibiting 
small anisotropy in the easy plane. They are described 
by the spin-1/2 XYZ spin model. By the appropriate 
design of molecules, it is possible to obtain a variety of 
spin systems, including those with spin larger than 1/2. 

The thermodynamics of the XYZ model has been 
solved exactly for the spin-1/2 case without external 



magnetic fields; in this case, the model is integrable and 
its HFE is obtained from the Bethe ansatz through a 
couple of integral equations^ and an infinite non-linear 
coupled equation[ll|. 

The model is no longer integrable, though, if an exter- 
nal magnetic field in any direction is introduced. This is 
the case of the XXZ spin-1/2 spin chain with transverse 
external magnetic field; so far, only its phase diagram at 
zero temperature has been investigated^, since the ex- 
ternal magnetic field does not commute with the Hamil- 
tonian. However, an analytic high temperature series of 
its HFE has not been obtained so far. 

In this communication we present the high tempera- 
ture thermodynamics of the spin-S" XY Z model, for ar- 
bitrary value of S, in the presence of an external magnetic 
field and a single-ion anisotropy term. We also investi- 
gate if the anisotropy in the XY plane avoids or does 
not avoid the classical behavior of the magnetic suscep- 
tibility and magnetization in this region of temperature, 
even for S = 3/2, as was shown in Ref. || for the spin-S 1 
XXZ model. 

The Hamiltonian of the spin-S XYZ quantum periodic 
chain with N sites reads 



iV 



H=£[J^SfS? 



ji qV qV 



jl qz qz 



h'Sf + D'iSlf 



(1) 



where Sf, a S {x,y, z}, are the spin matrices at the 
i-th site, the J' a are the exchange interaction couplings 
between first neighbors, and hi is the external magnetic 
field along the z-axis. We also include the single ion- 
anisotropy D' parallel to the external magnetic field. The 
periodic boundary condition S% +1 = Sf is used. 
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In order to render the thermodynamic functions to be 
finite, even in the classical limit (5 — > oo), we define 
a scaled spin operator s with unitary norm[s| as s e 
S/a/5(5 + 1). In order to write explicitly the effect of 
the anisotropy in the x and y directions in relation to the 
XX Z model 9], we rewrite the hamiltonian in terms 
of the spin operators sf and = -^(sj ± * s j); with 

N 

H = E{ J Mm+^ m +^ + + i 

i=l 

+s7 s - +l) + A s *s! +1 ]-h s ! + D( s l) 2 }. (2) 

The relations among the constants in hamiltonians Q 
and © are: J = (5/2)(5 + 1)(4 + J£), <5 ee - 

4)/( J * + J v)> A = 2 J 'J{ JL + 4)' h = y/W+T)h' and 



JJ ee 5(5 + 1)D'. 

In Ref. ^3 w e presented a closed version of the cum- 
mulant expansion for any chain (any one-dimensional 
classical or quantum model) with periodic boundary con- 
dition and interaction between first neighbours. A sur- 
vey with the main results of 0] is presented in Refs. 
|l4l . In this work we apply that method to compute 
the high temperature expansion of the HFE of the XY Z 
model (W s ), for a fixed value of the spin (with unitary 
norm). Like in the case of the XX Z model, this high 
temperature expansion is written as a series in powers of 
[5(5+ We use the interpolation method described 

in Ref. [9j to calculate the high temperature expansion 
of W s , for arbitrary values of 5, J, 5, A, h and D, up 
to order (J/3) 5 . The whole expression is too large, so we 
present it here only up to order (J/3) 3 , 
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(3) 



where h = h/J and D = D/J. 

We point out that these coefficients of the HFE, cal- 
culated up to order (J/3) 5 , are exact and valid for 5 = 
1/2, 1, 3/2, 2, • ■ • . The dependence of the HFE on even 
powers of S comes from the symmetry on the x and y 
directions. Letting Ws be the HFE of the XYZ model 
of spin with norm 5(5 + 1), we have that 

W S (J, 6, A, h, D- /3) = W s (ctJ 6, A, yfdh, aD; 0) , (4) 

where a = [5(5 + 1)] _1 . The HFE is an homogeneous 
function of first degree, so it is simple to obtain the /3- 
expansion of Ws(J, S, A, h, D; 3) from eq. ©• For 5 = 
we recover the results of Ref. 9] . The full expression of 
the expansion, up to order (J/3) 5 , can be obtained upon 



request to the authors. 

It is simple to derive from (j3J) the high temperature ex- 
pansion of the specific heat per site of the spin-5 of the 
XYZ model, with unitary norm (C s = -/3 2 ^f^ )- 
We obtain C s = - f + - f - | - 

2 |-)/3 2 + C(/3 3 ), which shows that in the high tempera- 
ture region, the XYZ model also presents a tail of the 
Schottky peak 16 (C'sch oc /3 2 ), for all values of 5. 

As a check of our /3-expansion of the HFE , valid for 
arbitrary spin-5, we compare the specific heat per site for 
S=l/2 derived from it to the numerical result obtained 
from the coupled equations by Takahashi[llT|. In Fig^we 
plot the specific heat for 5=1/2 with J x = — 1.0, J y = 
1.2 and J z = 2.0 (which corresponds to J = (3/4) x 0.1, 
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A = 20 and 5 = — 11) in the absence of an external 
magnetic field (h = 0) and no single-ion anisotropy term 
(D = 0). At T = 0.90 the difference of the solution is 
2.9%. 
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FIG. 1: Specif heat of spin-1/2 XYZ chain with h = 0, 
D = 0, J x = -1.0, J y = 1.2 and J z = 2.0. The solid line 
corresponds to the numerical solution of Takahashi and the 
dashed line represents the /3-expansion of this function derived 
from 101 . 



One difficulty about the spin-S* XYZ chain model is 
that it is no longer exactly soluble, in the presence of an 
external magnetic field, even for 5* = 1/2; moreover, the 
absence of symmetry in the x,y and z directions makes 
numerical calculations much more involved. However, its 
classical limit, in the high temperature region, can be 
easily obtained from eq.® by taking S — ► oo. 

In Fig[21we show that the magnetic susceptibility per 
site of this model can be approximated by their classi- 
cal behavior even for S — 3/2, in the high temperature 
region, with a percental error larger than 2%. 

FiglHI shows the comparision of the classical magneti- 
zation per site of the model to its quantum version for 
several values of spin, in the region of high temperatures. 
We verify that this thermodynamical function can be well 
approximated, in the high temperature region, by its clas- 
sical result up to S — 1; the percental error, in this case, 
is smaller than 3% (see Fig. 3.2). 

In summary, we have presented the high temperature 
expansion of the HFE of the XYZ model, for arbitrary 
values of spin, in the presence of an external magnetic 
field and single-ion anisotropy term, up to order (J/3) 5 . 
The thermodynamic functions derived from eq. © can be 
used to fit experimental data to determine the value of 
the constants that describe the material under interest. 
As a check, we show that our expansion of the specific 
heat coincides with the numerical solution of Takahashi's 
coupled equations [TH up to T ~ 1 for h = 0, D = 0, 
J x = -1.0, J y = 1.2 and J z = 2.0. 




FIG. 2: In 2.1 we compare the classical (S — + oo) and the 
quantum (S = 1/2, 1,3/2 and 4) magnetic susceptibility per 
site as a function of (J/3). In 2.2 we present the relative per- 
cental difference between the quantum and classical results. 
We let A = 1, 5 = 1, h/J = 0.3 and D/J = -0.5. 
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FIG. 3: The magnetization versus ft/J at J/3 = 0.4 is pre- 
sented in Fig. 3.1 for S = 1/2, 1,3/2,4 and S -* oo. Figure 
3.2 shows the percental difference of those quantum magneti- 
zation curves with respect to the classical one. We let A = 1, 
5 = 1, h/J = 0.3 and D/J = -0.5. 



We easily obtain the /3-expansion of the classical be- 
havior of the model in this region of temperature. Finally 
we showed that the magnetic susceptibility and magneti- 
zation of the quantum XYZ model can be approximated 
by their classical analog for S > 3/2 and S > 1 respec- 
tively. Our result allows the determination of the relative 
percental error between classical and quantum solutions, 
for any value of spin, for those thermodynamical func- 
tions, although this is not true for other functions like 
the specific heat per site. 
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